Simplified Topological Invariants for Interacting Insulators 
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We propose general topological order parameters for interacting insulators in terms of the Green's function at 
zero frequency. They provide an unified description of various interacting topological insulators including the 
quantum anomalous Hall insulators and the time reversal invariant insulators in four, three and two dimensions. 
Since only Green's function at zero frequency is used, these topological order parameters can be evaluated 
efficiently by most numerical and analytical algorithms for strongly interacting systems. 
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Introduction. Topological insulators are new quantum 
states of matter whose characteristic property is the exis- 
tence of both bulk energy gap and stable surface states] 1- 
@]. The stability of surface states are protected by the bulk 
topology, which in the non-interacting limit is described 
by Bloch band topological invariants such as the Thouless- 
Kohmoto- Nightingale-den Nijs(TKNN) invariant|5] and the 
Z2 invariants [ 6- 11]. More recently, topological insulator 
with strong electron-electron interaction is becoming a cen- 
tral topic in the field lU2H22ll . For general interacting sys- 
tems, the topological order parameters can be defined as the 
physical response function for the quantum Hall effect ll23ll 
and the topological magneto-electric effect] 11]]. For actual 
evaluations of these physical response functions, it was pro- 
posed that Green's function is an useful tool in topological 
insulators [24], and there is much recent interest focused in 
this direction|25-28]. However, our original formula for the 
topological order parameter !124l is rather complicated; more 
recently, a much simpler formula was obtained for the inver- 
sion symmetric interacting topological insulators] 29 1. 

The main purpose of this paper is to obtain several simple 
and yet general topological order parameters for interacting 
topological insulators in an unified framework. They are ex- 
pressed in terms of Green's function at zero frequency instead 
of the entire frequency domain. These invariants strongly 
resemble the conventional topological invariants such as the 
Chern number/TKNN invariant; yet, they are valid for gen- 
eral interacting systems. Current proposals for the quantum 
anomalous Hall (QAH) insulators| 30-32] require magnetic 
order, which is only possible for interacting systems. Our pro- 
posed topological order parameter can greatly help the search 
for realistic materials. Among our central results are Eq.(j6|, 
Eq.tO, Eq.dT6b, Eq.(Q~8]i, and Eq.(fT9b, all of which are ex- 
pressed in terms of Green's function at ico - 0. In most nu- 
merical algorithms of strongly interacting systems, it is much 
easier to obtain the Green's function at zero frequency rather 
than all frequencies. Therefore, our new formulas present 
a significant improvement over the previous result|]24|. We 
would also like to point out that the formulas given in this 
paper is not directly applicable to fractional topological insu- 
lators with nontrivial ground states degeneracy, which will be 
left to future studies. 

Topological order parameter for interacting QAH insu- 



lators. The conventional topological invariant for two- 
dimensional(2d) non-interacting quantum (anomalous) Hall 
states (or "Chern insulator") is the TKNN invariant[5], which 
is also called the first Chern number in mathematic literatue. 
Explicitly, it is an integral over the momentum space (namely 
the first Brillouin zone) 
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where fy - didj—djOi, and a,- = —i YiaW(^)\^K \^"(k)), where 
a runs through all the occupied bands. However, because its 
fundamentally dependence on the Bloch state \ip"{k)), it ap- 
plies only to non-interacting systems. There is an interesting 
generalization to interacting systems using twisted boundary 
condition] 33], which is nonetheless difficult to compute and 
is not easy to generalize to Z2 insulators. There is another in- 
teger topological invariant expressed in term of Green's func- 
tion rather than the Bloch state]24, 34, 35] 

N 2 = — - fdk^kTile^PGd^GdyG^GdpG- 1 ] (2) 

where fj., v,p run through k^, k\,k2, with ko = ico referring to 
the Matsubara frequency (imaginary frequency). Throughout 
this paper, the Green's functions are the Matsubara Green's 
functions, though our final formulas are also applicable using 
real frequency. For instance, a non-interacting system with 
Hamiltonian H - 2fcctA(&)c* has G(ico,k) - 1 /(ico - h(k)), 
where h(k) is generally a ^-dependent matrix, and q is a col- 
umn vector of fermion operators. The discrete Matsubara 
frequency becomes continuous in the zero temperature limit 
which we take. 

Eq.(j2]l has the severe disadvantage that it involves a fre- 
quency integral. In most numerical algorithms, it is very diffi- 
cult to obtain the dynamic Green's function at all frequencies 
as required by Eq.(f2}. Now we shall show that it is possible, 
without any approximations, to evaluate Eq.© with only the 
Green's function at zero frequency. Our new formula is much 
easier for practical calculations, and is accessible for most nu- 
merical algorithms. 

Let us start from the formalism presented in our previous 
work|29|. We diagonalize the inverse Green's function as 



G (ia>, k)\a(ioj, k)) — p. a (iu>, k)\a(iu>, k)) 



(3) 



The eigenvectors of G are the same as those of G l , with 
eigenvalues ji~ , therefore, we can also formulate our ap- 
proach by diagonalizing G instead of G ' . From the Lehmann 
representation, we can show that Green's function satisfies the 
equation 

(G- l )\iu>,k) = G-\-iu>,k) 



from which it follows that 

(G- l ) f (0,k) = G- l (0,k) 



(4) 



(5) 



Therefore, fi a (Q,k) are real numbers. The eigenvectors 
|ar(0, k)) can be divided into two types: those with /J, a (0, k) > 
are called "Right-zero (R-zero)", while those with /v„(0, k) < 
are called "Left-zero (L-zero)". All the R-zeros span a sub- 
space at each k, which we call the "R-space". Similarly we 
can define the "L-space". Because the eigenvectors corre- 
sponding to different eigenvalues of a Hermitian matrix are 
orthogonal, the R-space has the crucial property that all the 
vectors within it are orthogonal to those within the L-space. 
Therefore, the Chern number/TKNN number can be defined 
for the R-space in the conventional way. More explicitly, we 
present one of the central result of this paper, which we shall 
call the "generalized TKNN invariant", or the "generalized 
Chern number" 
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where T u = d t Jlj - djM t , andJfl< = -iT ia ER- sp ace(ka\d kt \ka). 
Here \ka) is an orthonormal basis of the R-space. The sim- 
plest basis choice is \ka) - \a{iw - 0, £))(i.e. normalized 
R-zero). Throughout this paper, the lower case expressions 
c\,C2,di,fij are reserved for the non-interacting Chem num- 
bers, Berry connection and curvature, while the upper case ex- 
pressions C\ , C2, Mu Tij refer to the generalized Chern num- 
bers, generalized Berry connection and curvature for gener- 
ally interacting systems. 

The expression in Eq.© reduces to the conventional TKNN 
invariant in the non-interacting limit. However, it is essen- 
tially different in that it is defined in terms of R-zeros and R- 
space of the Green's function of a many-body system rather 
than the Bloch states of a non-interacting Hamiltonian. The 
definition of R-zero and R-space is valid in the presence of 
general electron-electron interaction, and thus this topologi- 
cal order parameter is a highly nontrivial generalization of the 
TKNN invariant to the interacting insulators, maintaining the 
elegant form of Chern number. Mathematically, Eq.© is the 
first Chern number of an U(N) fiber bundle, whose bundle at 
each k is exactly the R-space. Therefore, this quantity is a 
topologically invariant. 

One can also define an "effective Hamiltonian" h e ff(k) = 
— G -1 (0, k), and define R-zeros as eigenvectors of this "effec- 
tive Hamiltonian". This is just another language to formulate 
the same result, which we shall use later in the derivation of 
our formula. 

The characteristic physical observable of the 2d correlated 
QAH insulators is the Hall conductance 
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which can be measured in experiments. As we shall show in 
the next section, the coefficient C\ here is exactly the one in 
Eq.©. The simple formula Eq.© can be easily applied to 
the interacting quantum (anomalous) Hall systems with inte- 
ger quantum Hall effect. However, it should be mentioned ll24ll 
that this description cannot be directly applied to the fractional 
quantum Hall states which have nontrivial ground state de- 
generacy. The same limitation is also true for the frequency 
integral formula in Eq.©. We will leave the possibility of ex- 
tending our approach to cases with ground state degeneracy 
for future studies. 

Derivation of the formula. Eq.© is itself a topological 
invariant that can be applied to interacting insulators in 2d 
that break time reversal symmetry. In this section, we would 
like to show that it is indeed the quantum Hall conductance, 
namely, N2 = C\, in other words, Eq.© can be derived from 
Eq.©. Let us begin with the Lehmann representation of Mat- 
subara Green's function (in the zero temperature limit) 



.*)=D 



(0\c ka \m)(m\cU0) (m\ Cka \0)(0\chm) 



ioj - (E m - E ) 



W 

iio + (E m - E ) 



(8) 



where \m) are the exact eigenvectors of K — H - fiN (H is the 
many-body Hamiltonian, and fi and N is the chemical poten- 
tial and fermion number respectively) with eigenvalues E m , 
and |0) is the ground state. We can make the summation over 
\m) well defined by putting the system on a large but finite two 
dimensional torus so that the eigenvalues are discrete, other- 
wise the summation is replaced by integral. Note that we have 
assumed that there is only a single ground state. For our pur- 
pose, we decompose G(ia>, k) — G\ + iGi with both Gi and G2 
Hermitian. Explicitly, we have 



(G2)a0 ~ ~ / , 



a ) 2 +(E m -E ) 2 
X[<0|cto|m><m|cj,|0> + <m| Cto |0><0|c^|m)] 
= -2_jd m [u ma {k)u ml3 {k) + v* ma (k)v m/} (k)] (9) 

m 

where we have defined u ma (k) = (m|cjj0), v ma = (0|cjjm) 
and d m (i(jL>) - (l>/[u> 2 + (E m - £o) 2 ] to simply the expressions. 
It is easy to see that sign(<i m ) = sign(w). Now we can calculate 
the expectation of G2 with an arbitrary vector \a) as (a\G2\a) = 

Tjap a*a(G2)al3 

this we can see 



" - Zm d m \\ Za a a u ma \ 2 +\ Z a a a v ma \ 2 ]. From 



sign((a|G2(Jo>, k)\a)) - -sign(w) 



(10) 



As a side remark, if \a) is an eigenvector of G = G\ + /G2 with 
eigenvalue /j~ l , then we have /i^ 1 = (a|a) _I (fl|(Gi + iG2)\a). 
Due to the fact that G\ and G2 are Hermitian, we have 
lm(p. a l ) - (a|fl) _I (fl|G2|a), thus it follows from Eq.([T0t that 
sign[lm(ju~ (iaj))] - -sign(w). It follows finally that 



sign[lm(ju a (/w, k))] = sign(w) 



(11) 



(7) 



With these preparations, we are approaching the central part 
of our calculation. The key idea is to introduce a smooth de- 
formation of G(icu), k) parameterized by A e [0, 1] as follows 

G(ico, k,A) = (1 - A)G(iaj, k) + A[icj + G l (0, k)y l (12) 



We now show that this deformation does not contain singu- 
larity, or equivalently, we must show that all eigenvalues of 
G(iu>, k, A) are nonzero. This can be seen as follows. First, 
when ico - 0, we have G(0, k, A) - G(0, k), whose eigenval- 
ues are nonzero by our assumption that the Green's function 
G(ia>, k) is nonsingular. Second, we consider ico + 0. Suppose 
that 

Giico, k, A)\a(ico, k, A)) = li~ {ico, k, A)\a(ico, k, A)) 
then we have 

^(ico, k, A) - {a\ay l {a\G(iu>, k, A)\a) 
The imaginary part of this equation can be written down as 

Im[p a l (ico, k, A)] = (a\a)- l [(l - A)(a\G 2 (ico,k)\ a ) 

-aoj^sHoj 2 + 4r l ] 

s 

where we have expanded \a(ico, k, A)) - J^a'sOw, k, A)\s(k)}, 
in which \s{k)} are orthonormal eigenvectors of -G~'(0, k) 
with eigenvalues e s (k). It is easy to see that lm\ji~ l (iu>, k, A)] 
is always nonzero thanks to Eq.dTOb. Summarizing the above 
calculation, we can see that all eigenvalues of G(ia>, k, A) are 
nonzero, therefore the deformation in Eq. dTZi is smooth. Note 
that throughout this calculation we consider the imaginary fre- 
quency Green's function, otherwise Green's function cannot 
be so well-behaved. A geometrical visualization of the de- 
formation in Eq.dT2b can be given as follows. Because of 
Eq.dTTI). the p a (ico)(co e (-00, +00)) curves[29] on the com- 
plex plane do not cross the real axis when ico + 0, therefore 
we can smoothly deform them to straight lines parallel with 
the imaginary axis, keeping the R/L-zero unchanged in the 
deformation. This leads exactly to Giico, k,A- 1). 

Because N 2 is a topological invariant, namely that it is un- 
changed under smooth deformations of G, we have N 2 (A = 
0) = N 2 (A = 1). Therefore, to calculate N 2 = N 2 (A = 0), 
we just need to calculate N 2 (A - 1), which is equivalent 
to the calculation for an effective noninteracting system with 
h e ff(k) = -G~'(0, k). It is a straightforward calculation to 
obtain N 2 (A - \) - C\. This completes the derivation of 



N 2 = C X 



(13) 



which is a precise identity between Eq.(f2]i and Eq.©. 

Four-dimensional topological insulators. The 2d physics 
discussed above can be generalized to 4d. In 4d, there is 
time reversal invariant topological insulator classified by in- 
teger Z. The continuous model for such topological insulators 
was first proposed in Ref.[36], while lattices models can be 
found in Ref. lUlll . For non-interacting insulators in 4d, the 
natural topological invariant is the second Chern number in 
the momentum space expressed as[ 1 1] 



with 
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L- 2 J d 4 ke^tr[f u f kI } (14) 
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af(k) 



n a/3 a ail , .f V 

oia. - Ojdj + 1 ai, tZj 



where i, i,k,l — 1,2,3,4. The index a in a- refers to the 
occupied bands of the Bloch states \if/ a (k)}. The Berry con- 
nection a" is a non-abelian gauge field potential, and fr is 
the associated non-abelian field strength. Analogous to the 2d 
case, for 4d interacting insulators there is a topological invari- 
ant analogous to Eq.(f2|, expressed in terms of the interacting 
Green's function II 111,12411 

^4 r yj^ f d^Trle^f^Gd^GdyG^GdpG- 1 

A Gd T G- 1 ] (15) 



480tt 3 
xGdM- 



This topological order parameter directly measures the gen- 
eralized quantum Hall effect in 4d llll |36||. and it is related 
to the homotopy group 7T5(GL(A^, C)) — Z[24] for sufficiently 
large N. The difficulty with Eq.dT5b is again the frequency 
integral over (— ioo, +ioo). This problem can be solved in the 
same way as its 2d analog. We are thus led to another central 
result of this paper, namely the topological order parameter 
for an interacting Chern insulator in 4d, expressed as 



(16) 



with 



ii/j 



t* = dirt* - dj&f + i [& u &j] 

Mf(k) = -i(ka\—\kj3) 

okj 

where \ka) is the same as that defined below Eq.©, namely 
that I to) is an orthonormal basis of the R-space spanned by 
R-zeros. The derivation of N4 — C 2 is a straightforward gen- 
eralization of that of Eq.([T3l. which we shall not repeat here. 
Z2 topological invariants for interacting insulators in three 
and two spatial dimensions. It was proposed in Ref.[24] that 
a natural topological invariant for 3d X 2 insulator is the topo- 
logical magneto-electric coefficient ll24l 12911 . 
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xGdyG^GdpG^Gdo-G^GdrG 1 } 



(17) 



where the integer W(G)\ RxT 4 is the "winding number" of the 
mapping from frequency-momenta space R x T A to GL(N, C), 
ko = ico is the imaginary frequency, and k\ is the dimensional 
extension parameter similar to the Wess-Zumino-Witten pa- 
rameter in non-linear cr models. The reference function 
G(ko,ki,k 2 ,k2,7i) is trivially diagonal ll24ll . Due to the am- 
biguity of the dimensional extension, the integer in Eq.dTTb 
reduces to Z 2 equivalent classes[24]. For insulators with in- 
version symmetry, Eq.([P71 is further simplified to a product of 
parity of R-zeros [29]. This major simplification enables prac- 
tical numerical calculations using Green's function, see e.g. 
& 

Eq.dTTb. although an elegant Z 2 invariant, is unsatisfactory 
because the needs of the dimensional extension and the fre- 
quency integral. Now we shall obtain Z 2 topological invari- 
ants without these two disadvantages. This is outlined as fol- 
lows. Let us start from Eq.dTTb and dimensionally extend 



G(ko,ki,k2,kj,) to G(ko,ki,k2,k3,k4), where k\ e [—n,ri\ is 
the WZW like dimensional extension parameter. Zero fre- 
quency function G(0, k\,ki, £3, £4) are chosen to be Hermitian, 
therefore, R-zeros can be extended to the extended momenta 
space T 4 . By a calculation analogous to the the derivation of 
N4 = C2, Eq.dT7b can be simplified into an interacting topo- 
logical order parameter in 3d as 

P3 = C 2 /2 =-L^fdk 4 f dh^'tTlTuTk,] 

= JL f dh^Triid^jik) + -mmftj-mmkik)} us) 

Here the Berry connection is defined in terms of zero fre- 
quency Green's function in the same way as Eq.© and 
Eq.(fT6b. Eq.(fT8l) generalizes the formula first obtained by Qi, 
Hughes and Zhang for the non-interacting system lTllll . Be- 
cause the Green's function instead of the Bloch states are used, 
Eq.dTSb is defined for interacting topological insulators in 3d. 
Eq.dTSb has great advantage compared to Eq.(fl7t because nei- 
ther dimensional extension nor frequency integral is needed. 
In the non-interacting limit, the R-zeros become Bloch states, 
and Eq.(fT8l) is thus reduced to the Chern-Simons(CS) invari- 
ant defined in Ref. lUlll . As was shown in Ref.[38], the non- 
interacting CS invariant! 11] is equivalent to the Fu-Kane's 
Pfafhan invariant] 8], it is thus natural to obtain an Pfaffian 
invariant A for interacting topological insulators, which is ex- 
pressed in terms of R-zeros as 



(-d a = n 



r /= TRIM 



VdetB(r,) 
PW/)) 



(19) 



Here, "TRIM" refers to time reversal invariant momenta, the 
matrix B(k) is defined by B a p - (-ka\T\kB), where T is the 
time reversal operation, and \ka) is an orthonormal basis of 



R-space (the simplest choice is \ka) - \a(iu> = 0, k))). Eq.([T8b 
and Eq.(fT9t, which are Z2 invariants defined for interacting in- 
sulators, are among the central results of this paper. They have 
the great advantage that only zero frequency Green's function 
is needed. The ambiguity of the square root in Eq.(|T9b can 
be avoided if one rewrite it as an integral following Ref.[8J. 
We also mention that there are other formulas equivalent to 
Eq.([T9l, e.g. number of zeros of Pfaffian in the half Brillouin 
zone, which parallel the non-interacting formulas, with Bloch 
states replaced by vectors in R-space. The derivation from 
Eq.(fT8]l to Eq.(fT9l> follows exactly the calculations in Ref.[38], 
again with the Bloch states replaced by vectors in the R-space. 
A similar expression as Eq.(|T9b can also be defined for 2d in- 
teracting topological insulators. For insulators with inversion 
symmetry, starting from Eq.([T8~t or Eq.dT9l one can derive the 
parity formula in Ref.[29], which was originally derived from 
Eq.dTTb directly. The parity formula 112911 is most convenient 
for practical calculation if the insulator has inversion symme- 
try. However, if the insulator has no inversion symmetry, we 
need to use the more general formulas in Eq.([T8l or Eq.([T9} 
proposed in the present paper. 

Conclusion. In this work we present a general framework to 
describe interacting insulators in terms of the Green's function 
at zero frequency. Our central results include the topological 
order parameters for the quantum anomalous Hall insulator, 
interacting Landau level systems with integer quantum Hall 
effect, and time reversal invariant interacting topological in- 
sulators in 4d, 3d and 2d. These formulas greatly simplify 
numerical and analytical calculations. 
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